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Abstract 

In a recent paper Goriely considers the one-dimensional scalar reaction- 
diffusion equation m = u xx + f(u) with a polynomial reaction term f(u) 
and conjectures the existence of a relation between a global resonance of the 
hamiltonian system u xx -\-f{u) = and the asymptotic speed of propagation of 
fronts of the reaction diffusion equation. Based on this conjecture an explicit 
expression for the speed of the front is given. We give a counterexample to 
this conjecture and conclude that additional restrictions should be placed on 

the reaction terms for which it may hold. 
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I. INTRODUCTION 



The one dimensional scalar react ion- diffusion equation 

u t = u xx + f(u) (1) 

with /(0) = f(u + ) = has been the subject of much study not only because it models 
different phenomena [|T]-0], but also because it is the simplest reaction-diffusion equation 
for which rigorous results can be obtained . Depending on the situation being 

considered the reaction term f(u) satisfies additional properties. It has been shown for 
different classes of reaction terms that suitable initial conditions u(x, 0) evolve in time into 
a monotonic front joining the state u = u + to u = 0. The asymptotic speed at which the front 
travels is the minimal speed for which a traveling monotonic front u(z) = u(x — ct) exists 
f5|j7|]. Traveling fronts are a solution of the ordinary differential equation u zz + cu z + f(u) = 0. 
In the present case we shall be concerned with two types of reaction terms, the classical case 
/ > in (0, u+) with /'(0) > and the bistable case / < in (0, a), / > in (a,u+) with 
Jq + f > and /'(0) < 0. In the classical case there is a continuum of speeds c > c* for 
which monotonic fronts exist. The system evolves into the front of speed c*. In the bistable 
case there is a unique isolated value of the speed c* for which a monotonic front exists, the 
system evolves into this front. The problem is to determine the speed of propagation of 
the front. In the classical case, if in addition /'(0) > f(u)/u the speed of propagation c* is 



the so called linear or KPP value ckpp = 2y/' (0) §]. In the other cases (as well as in the 
bistable case) there exist variational principles both local and integral from which the speed 
can be calculated with any desired accuracy for arbitrary / ]3|,| , H| lj] . 



In a recent paper |12| Goriely proposes a new method for the determination of the speed. 
Based on an observed property of some exactly solvable cases, namely reaction terms of the 
form f(u) = /iit + vu n — w 2n-1 he conjectures that for polynomial reaction terms of the form 

f(u)=nu + g(u) (2) 



2 



where g(0) = g'(0) = and the polynomial g independent of //, the speed of the front can 
be calculated from the knowledge of the heteroclinic orbit of the Hamiltonian system 

u zz + f(u) = 0. 

This property of solvable cases had not been observed before. The purpose of this article 
is to show by means of a counterexample that this conjecture is not true in general for 
the above class of polynomial reaction terms. However, considering that the conjecture is 
indeed true for a large class of reaction terms (the solvable cases mentioned above), it is 
a very interesting unsolved problem to characterize the class of functions for which it is 
valid. For the sake of clarity we state the conjecture here. The conjecture makes use of the 
fact the front approaches the equilibrium state u = as e A ~ 2 , a well established fact, and 
approaches the equilibrium point u = u + as u = u + — Le 7+Z , an assumption which is not 
always satisfied. Then the global resonance, defined as 

(3) 

is conjectured to be a constant, for a general class of polynomial reaction terms, at all values 
of fj, for which the nonlinear front exists. Explicit expressions are known for the rates of 
approach 7+ and A_ in terms of c and / therefore, if 5 can be calculated at any point, 
then an analytic formula for the speed can be obtained. There is such a point where it 
can be calculated, and that is the point at which c = and the system is hamiltonian. 
There is a unique value /i < for which such a front exists, (we shall label it as /i^ and the 
corresponding equilibrium point as Uh) and therefore the speed is completely determined. 
In the following section we consider a specific polynomial reaction term and show that it 
fails to satisfy the conjecture. 



II. THE COUNTEREXAMPLE 

Consider the reaction term f(u) = fiu + 2u 2 — 7u 3 + y-u 4 — 2u 5 . This is of the form 
f(u) = [ill + g(u) where the polynomial g(u) = 2u 2 — 7u 3 + 20w 4 /3 — 2u 5 satisfies the 

3 



properties g(0) = g'(0) = and is independent of /x as requested by the conjecture. For the 
hamiltonian system 

20 

u zz + fm+ 2u 2 - 7u 3 H u 4 - 2u 5 = 

3 

a heteroclinic orbit joining two equilibrium points exists at the value /x — Hh = —0.153897. 
In Fig. 1 the reaction term f(u) is shown together with the (scaled) potential. It is clear 
that a heteroclinic solution joining the point u = to u + = Uh = 0.262156 exists. In this 
case the resonance 5 can be calculated. Its value is given by 



5 = 5 h = Ji^- = 0.865558. (4) 

V 

Let us now consider the propagating fronts which are a solution of 
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Before giving the results of the numerical and analytical calculations we show the plot of the 
function / at several values of /x which will make clear the numerical and analytical results 
that follow. As /x increases the equilibrium point u + increases until at /x = 1/3 it reaches 
the value u + = 1 where /' = 0. Above this value of /x there is a discontinuous jump in u + , 
the front joins the origin u = to a new fixed point which corresponds to a different root of 
the polynomial f{u). In Fig. 2 we show the function / at different values of fi. At /x = 1/3 
the fixed point u + = 1 and the derivative f'(u + ) = 0. At /x = 0.4, we see that the value of 
u + is now the new root of / which did not exist at low values of /x. 

First we describe the results of the numerical integrations of the initial value problem for 
Eq(|I|) with sufficiently localized initial value perturbations u{x, 0). The speed is obtained 
numerically and the value of S is then computed from Eq. (R|) which can be expressed as W2 



-c + Jc 2 - 4/'(u+) 

6 = V ; 1 . (5) 

c+ / C 2_ 4f(0 ) 



In Fig. 3 we show the asymptotic speed of the front as a function of /x. The solid line gives 
the numerical results and the dashed line corresponds to the linear or KPP value 2-JJi. First 
we observe that the KPP value is lower than the calculated speed in the range of /x shown 



which means that the transition to the linear or KPP regime occurs at larger values of p. 
Even though the value of u + is discontinuous, the speed is a continuous function of p. 

And finally the graph of the resonance 5 as a function of p is shown in Fig. 4. At p = ph 
it adopts the analytically calculated value from the hamiltonian case, decreases to a value 
<5 = Oat/i = l/3 jumps discontinuously to a larger value and increases from there on. This 
discontinuity can be attributed to the discontinuity in u+. At the value of p = 1/3 where 
f'{u + ) = it is evident from Eq (|5|) that 5 = 0. As we will show below, at this value of p 
the speed and the asymptotic behavior for the front can be calculated analytically and it is 
found that the front does not approach the fixed point u + = 1 exponentially, therefore one 
of the assumptions of the conjecture does not hold. Indeed we shall show that at /i = 1/3 
the front approaches u = 1 as«~ 1 — A/z. 

In order to determine the speed we shall make use of variational principles. It is known 
that the speed of the front is given by 



c = max 2 (6) 

J + (j)du 

where the maximum is taken over positive decreasing functions <f)(u). Taking as a trial 
function 

= (l-^) 7/2 e -3/[(2(l-u)] 
>U 



the integrals in Eq.(||) can be performed. We obtain for /x = 1/3 

c > ^ > 2^/1 

To obtain an upper bound we make use of the local variational principle 

c = inf sup { p'iu) + ^ . . ) 
p u \ p{u)J 

where the trial function p{u) > and p'(0) > 0. Choosing as a trial function p(u) 



2/3 u (1 — u) we obtain the upper bound 



which combined with the above lower bound implies that the speed is exactly c = y 3/2 in 
agreement with the results of the numerical integration. The exact value of the speed could 
be obtained analytically from the variational principles due to the fact that for \i — 1/3 the 
derivative of the front can be calculated exactly. The derivative of the front as a function of 
u, p{u) = —dujdz satisfies the equation p{u)p'{u) — cp{u) + f{u) = and the exact solution 
at fi = 1/3 is given by p(u) = J2/3u (1 — u) 2 . With this expression for p we may calculate 
the approach to the fixed point u—1. Near u — 1, p ~ v/2/3 (1 — u) 2 so that 




from where it follows that 



u(z) ~ 1 — 




We see then that at this point 5 = since the rate of approach is not exponential but 
algebraic and one of the assumptions of the conjecture is not satisfied. We conclude from 
this example that the conjecture does not hold for general polynomials of the form given by 
Eq.@. 



III. CONCLUSION 

We have seen by means of a counterexample that the conjecture put forward that relates 
certain properties of the hamiltonian system u zz + f(u) = with the speed of the front 
solution of u zz + cu z + f{u) = is not satisfied by general polynomial reaction terms f{u). 
As observed by Goriely, there is a class of reaction terms for which it does hold, those for 
which an exact solution for the front u(z) can be given explicitly. Numerical evidence has 



been given [|T^] in at least one case where the conjecture seems to hold in a case where the 
front cannot be calculated explicitly. On the other hand, we have given a counterexample 
to this conjecture. It is an interesting problem to establish precisely the conditions under 
which the proposed conjecture holds. This would lead to a classification of systems in at 
least three classes, those for which the speed is given in terms of the derivative at one fixed 
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point, that is the KPP value c = 2y/'(0), those in which the speed would be determined 
by the derivatives at the two fixed points (the expression given by Goriely would hold) and 
the rest, for which the speed depends on integral properties of the reaction term. 
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Figure Captions 



Figure 1 

Graph of the reaction term / and the corresponding scaled potential at the value of fi 
for which the speed of the front vanishes and the system is hamiltonian. 

Figure 2 

Graph of the reaction term at different values of \x. The value of the stable point increases 
with /i until ji reaches 1/3. A discontinuous jump in the stable point occurs at that value. 

Figure 3 

Graph of the speed obtained from the numerical integration of the initial value problem. 
The speed of the front is a continuous function of /i. In the range of /x shown the speed is 
greater than the linear or KPP value. 

Figure 4 

Value of the resonance 5 as a function of ji obtained from the numerical integrations. 
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